Charge order in Fe 2 OB0 3 : An LSDA+[7 study 
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Charge ordering in the low-temperature monoclinic structure of iron oxoborate (Fe2 0B03) is 
investigated using the local spin density approximation (LSDA)+(7 method. While the difference 
between ti g minority occupancies of Fe 2+ and Fe 3+ cations is large and gives direct evidence for 
charge ordering, the static "screening" is so effective that the total 3d charge separation is rather 
small. The occupied Fe 2+ and Fe 3+ cations are ordered alternately within the chain which is in- 
finite along the a-direction. The charge order obtained by LSDA+t/ is consistent with observed 
enlargement of the j3 angle. An analysis of the exchange interaction parameters demonstrates the 
predominance of the interribbon exchange interactions which determine the whole L-type ferrimag- 
netic spin structure. 



PACS numbers: 71.20.-b, 71.28.+d, 71.30.+h 

I. INTRODUCTION 

Transition metal compounds with interplay of spin, 
orbital, and charge degrees of freedom are of strong 
current interest. Thus, magnetite (FeaC^), the famous 
lodestoneji is one of the classical examples of such a 
system, where a first-order metal-insulator transition oc- 
curs at ~ 120 Ki£ According to Verwey this transition 
is caused by the ordering of Fe 2+ and Fe 3+ cations on 
the octahedral B-sublattice of the inverted spinel struc- 
ture AP^C^ii^ Recently, the local spin density approxi- 
mation (LSDA)+[/ method 5 '- 6 ^ has been used for inves- 
tigation of the low-temperature monoclinic structure 8 9 
of FeaO^iS. The charge and orbitally ordered insulating 
ground state with two different order parameters: small 
for the charge and large for the orbital order was obtained 
self-consistently. While the first order parameter was in- 
troduced as the total 3d charge separation between Fe 2+ 
and Fe 3+ cations, the latter one corresponds to the charge 
difference between occupancies of occupied t2 9 i orbital of 
Fe 2+ and unoccupied i2gi orbital of Fe 3+ cations. But an 
important issue, which is still unresolved is whether the 
charge ordering in magnetite is driven by Coulomb repul- 
sion between the charges, or by the strain arising from 
electron-lattice interactions. An analysis of the charge 
and orbital order parameters 1 ^ in Fea04 inevitably re- 
sults in a strong interplay between Jahn- Teller effect 
and electrostatic repulsion between electrons. Therefore, 
theoretical investigations of the isostructural iron and 
manganese oxoborate (Fe20B03 and Mn20B03, cor- 
respondingly) that support predominantly electrostatic 
and Jahn- Teller distortion driven mechanisms^ behind 
the charge ordering in the same structural arrangement, 
are highly desirable. 

Iron borate (Fe20BOa) is a semi-valent oxi&eM^ It 




FIG. 1: The Fe 2 OB0 3 structure projected on (100) plane (6 
vertical, c horizontal). Two structurally distinct Fe(l)06 and 
Fe(2)Os octahedra shown by light and dark sharing respec- 
tively. Plus and minus signs indicate the relative orientation 
of the moments within each Fe(l) and Fe(2) chains in the 
magnetically ordered phase. 



belongs to the homomctallic warwickite family with for- 
mal chemical formula MM'OBO^, where M and M' are, 
respectively, a divalent and trivalent metal ions. Surpris- 
ingly, that the homometallic (M — M') warwickites are 
known only for FeiSiii and Mnm^^ In both compounds 
the metal have octahedral coordination. These octahedra 
share edges to form ribbons of four infinite along crystal- 
lographic a-direction chains of octahedra linked by corner 
sharing and the trigonal BO3 groups (see Fig. QJ. 

There are two crystallographically inequivalent sites 
of the metal ions Fe(l) and Fe(2). Fe20B03 is L-type 
ferrimagnetic with drastically smaller in comparison to 
Fe3C>4 critical temperature of T c w 155 K, the Fe(l) 
magnetic moments being aligned antiparallel to the 
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Fe(2) moments. It is almost antiferromagnetic but a 
small fcrrimagnctic moment of ~ 0.03 /x^ per Fe atom in 
a 0.05 T field was found ,Iii±2ii&±£i At room temperature 
Fe20B03 is a semiconductor with a thermoactivated 
conductivity low e~ E °/ kT with E a ps 0.35 eViiiS Upon 
farther heating a broad semiconductor-to-semiconductor 
transition occurs at T co m 317 K, where resistivity 
drops down by a factor of ~ 3, and, as a result, a small 
decreasing of the activated energy up to E a m 0.31 eV 
above 350 K is observediiii^ The 317 K transition is 
assigned to charge ordering of 2+ and 3+ Fe cations on 
Fe(l) and Fe(2) sites, and accompanied by a structural 
transition from monoclinic P2\/c to orthorhombic 
Pmcn symmetry with increasing temperature. This 
structural transition is attributed by modification of 
the (3 angle from (3 = 90.220(1)° at 3 K to [3 = 90° 
at 337 Kiiiiiaii^ The change in conductivity and struc- 
ture are small. But the 57 Fe Mossbauer spectra at 
around 317 K clearly result in the charge localization 
at the transition with an equal distribution of Fe 2+ 
and Fe 3+ cations over the two structurally distinct 
Fe(l) and Fe(2) sites with formal chemical formula 
Fe(l)g + 5 Fe(l)g + 5 Fe(2)g+Fe(2)g+OB0 3 ai™saL^ 
Although, there are two types of distorted FeOg oc- 
tahcdra with Fe-0 bond length varying between 1.92 
and 2.23 A for 3 K, the average Fe(l)-0 and Fe(2)-0 
distances are 2.085 and 2.082 A, respectively, i.e. equal 
within experimental errors^ Such a small difference 
results in the extremely small value of deviation (< 0.01) 
from the average 2.5+ value of valence of Fe cations 
estimated by the bond valence sum method. While 
an electronic transition between charge ordered and 
disordered state occurs at around 317 K, as evidenced by 
Mossbauer spectroscopy and resistivity measurements, 
no long range Fe 2+ /Fe 3+ ordering is directly observed by 
x-ray, neutron or electron diffraction. Thus, a long range 
charge ordering such as the simple alternating scheme 
proposed in Ref. [llj destroys the mirror symmetry, 
which leads to a tilting of the Fe-ribbons, consistent 
with the observed enlargement of the (3 angle below 
the transition. However, there is no observation of the 
increasing of a-axis periodicity (it should increase by a 
factor of two or another integer factor below T co ). Thus, 
below the transition, a charge ordering is not implicit in 
the atom coordinates, although it is indirectly evidenced 
by other experiments. This ambiguity is resolved in our 
electronic structure study, which reveals an arrangement 
of Fe 2+ and Fe 3+ cations alternately ordered within the 
chains along the a-direction. 

In this paper we report theoretical investigation of the 
electronic structure and magnetic properties of Fe20BC>3 
in the low-temperature P2\jc structure. The LSDA+J7 
approach in the tight-binding linear muffin-tin orbital 
(TB-LMTO) calculation scheme^ has been used. Mo- 
tivated by our results, we propose an order parameter, 
defined as the difference between ti g minority spin oc- 
cupancies of Fe(l) 2+ and Fe(l) 3+ as well as the differ- 
ence between t^g majority spin occupancies of Fe(2) 2+ 
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FIG. 2: Total density of states (DOS) obtained from the 
LSD A calculations for the low-temperature P2\/c phase of 
Fe20B03. The Fermi level is shown by dotted line. 



and Fe(2) 3+ cations. This order parameter is found to 
be quite large, although the total 3d charge difference 
between 2+ and 3+ cations, is small. 



II. COMPUTATIONAL DETAILS 

The present band-structure calculations have been car- 
ried out for the low-temperature monoclinic structure of 
Fe20BC>3. The corresponding P2\/c unit cell contains 
four Fe20B03 formula units. We used atom coordi- 
nates from Ref. ^] (see Table 1) and the cell parame- 
ters a = 3.1688 A, b = 9.3835 A, c = 9.2503 A, and 
(3 = 90.22° refined at 3 K. n i 12 i 13 The radii of muffin- 
tin spheres were taken Rp c = 2.5 a.u., i?s = 1-275 a.u., 
-Roi,04 = 1-7 a.u., and i?02,03 = 1-8 a.u. Seven kinds 
of empty spheres were introduced to fill up the inter- 
atomic space. For simplicity we neglect small spin-orbit 
coupling (for instance, in other iron oxide, in Fe 3 04, the 
spin-orbital interaction for 3d electrons was found to be 
negligibly small) . We consider only a collinear spin case, 
in which the Fe magnetic moments are aligned along the 
a-direction (see Table 1 of Ref. Il3t ). which is in a rea- 
sonably good agreement with the experimental magnetic 
structure of Fe20B0 3 below T c . 



III. LSDA BAND STRUCTURE 

The LSDA calculations give only a metallic ferrimag- 
netic solution without charge separation where partially 
filled bands at the Fermi level originate from the t^g or- 
bitals of Fe cations (see Fig- • 

The lower part of the valence band (below -3.5 eV) is 
mainly formed by O 2p states with a bonding hybridiza- 
tion with Fe 3d states. Fe 3d states give predominant 
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contribution to the bands at -3.5 eV below and up to 
2.5 eV above the Fermi level. The exchange splitting be- 
tween the spin-up and spin-down Fe 3d states is roughly 
3 eV. Additionally, the five-fold 3d levels are split by the 
crystal field into t 2g and e g sub-bands. The oxygen oc- 
tahcdra in Fe2 0B03 are strongly distorted and the local 
symmetry of Fe sites is, of cause, lower than cubic. Nev- 
ertheless, the cubic component of the ligand field, which 
is determined by the relative strength of Fe d - O p hy- 
bridization of 7r- and (7-type, remains dominant, whereas 
the splitting within u t 2g " and "e g " subbands is smaller 
than the corresponding band width. This allows one to 
label the corresponding states as t 2g and e g . The crys- 
tal field splitting is roughly 2 eV, which is less than the 
exchange splitting. This is consistent with the high-spin 
state of the Fe cations. The symmetry inequivalence of 
Fe(l) and Fe(2) sites leads to an inexact cancellation of 
magnetic moments and results in a small ferrimagnetic 
moment of ~ 0.31 fj,B per formula unit. The absolute 
values of magnetic moments obtained by LSDA are 3.54 
fiB and 3.81 \j,b for Fe(l) and Fe(2) sites, respectively. 
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FIG. 3: The total DOS obtained from LSDA+17 calculations 
with 17=5.5 eV and J=l eV for the low-temperature P2i/c 
phase of Fe20B03. The top of the valence band is shown by 
dotted lines. 



Fe(l) and Fe(2) t 2g and e g states with the opposite 
spin projections share nearly the same energy intervals. 
Thus, Fe 3d states between -3.5 and -2.0 eV originate pre- 
dominantly from majority spin Fe(l) and minority spin 
Fe(2) t 2g states whereas the states between -2.0 and - 
0.5 eV are mainly of e g character. Partially occupied 
bands crossing the Fermi level are formed by minority 
spin Fe(l) and majority spin Fe(2) t 2g states. The nomi- 
nal occupation of these bands is 1/6. In the majority spin 
channel, however, the Fe(2) t 2g state, that is oriented in 
the plane perpendicular to the shortest Fe(2)-0 bond, 
forms quasi-one-dimensional bands with a strong disper- 
sion along the a-direction. The one-dimensional charac- 
ter of the dispersion is explained by the fact that there 
are only two nearest neighbours of the same kind around 
each Fe(2) ion. The other two Fe(2) t 2g states are shifted 
to higher energy and the corresponding bands arc com- 
pletely unoccupied. As a result, the majority spin bands 
crossing the Fermi level turn out to be half-filled. An 
Fe(l) ion, in contrast to Fe(2) one, has four Fe(l) neig- 
bours at close distances. As a result of the hybridisation 
between Fe(l) t 2g states the situation in the minority 
spin chancl is more complicated. Twelve t 2g bands are 
split into three groups of 4 bands each. The Fermi level 
is crossed by lowest bands which show a rather strong 
dispersion along a-direction but with a period two times 
smaller than the quasi-one-dimensional Fe(2) bands. 



It should be noted that in contrast to experimental 
data n i 17 i 18 LSDA predicts Fe 2 OB0 3 to be metallic with 
substantial magnetic moment per unit cell. Apparently, 
the electron-electron correlations, mainly in the 3d shell 
of Fe cations, play a significant role. 



IV. LSDA+(7 RESULTS AND CHARGE 
ORDERING 

To proceed further we take into account the strong 
electronic correlations in Fe 3d shell using the LSDA+J7 
method. The calculations have been performed for the 
P2\/c unit cell as well as for double (2a x b x c) and 
triple (3a x b x c) P2\/c supercells of Fe20BC>3 (without 
putting in any additional local displacements of oxygen 
atoms around of Fe 2+ /Fe 3+ sites). For the 2a x b x c su- 
percell the LSDA+J7 calculations using the classical value 
of Coulomb (Z7=5 eV) and exchange (J=l eV) interac- 
tion parameters for Fe result in an insulating charge or- 
dered (CO) solution with an energy gap of 0.13 eV. This 
is in a strong contrast to metallic solution without CO 
obtained by the LSDA. The CO pattern obtained from 
the calculations is similar to the one proposed in Ref. ^] 
This is a notable result because the CO is not implicit 
in the atomic coordinates, and it shows that LSDA+?7 
calculations can assist experiments in revealing CO ar- 
rangements. To obtain a reasonably good agreement of 
the calculated gap of 0.39 eV with experimental value of 
0.35 eV we increase the U value up to 5.5 eV (see Fig- - 
It does not exceed of 10% of the U value, which is in an 
accuracy of the U calculation. Note, however, that the 
CO obtained by LSDA+J7 within 2axbxc supercell does 
not depend on the U value of 5-5.5 eV. Here and in the 
following all results are presented for the double along 
a-direction P2i/c supercell of Fe 2 OB03. 

After self-consistency each of two groups of Fe(l) and 
Fe(2) atoms is split out in two subgroups of 2+ and 3+ Fe 
cations with equal number of 2+ and 3+ cations. Thus, 
one of t 2g majority/minority spin states of Fe(2)/Fe(l) 
atom becomes completely occupied, whereas all other t 2g 
states are pushed by strong Coulomb interaction at the 
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FIG. 4: The partial DOS for different Fe cations are shown. 
The gap is opened between Fe(2) 2+ and Fe(2) 3+ for majority 
spin and Fe(l) 2+ and Fe(l) 3+ cations for minority spin states. 
The gap value of 0.39 eV was obtained by LSDA+C7 with 
£7=5.5 eV and J=l eV. The Fermi level is shown by dotted 
line. 



energies above 3 eV. The gap is opened between occu- 
pied and unoccupied ti g states of Fe(l) 2+ and Fe(l) 3+ 
for spin-down and Fe(2) 2+ and Fe(2) 3+ for spin-up. Ma- 
jority spin 3d states of Fe(l) 3+ and minority spin states of 
Fe(2) 3+ cations are shifted below the O 2p states, which 
form the band in the energy range of -8 and -2 eV. In 
contrast to Fe 3+ states, the majority spin Fe(l) 2+ and 
minority spin Fe(2) 2+ 3d states form the broad bands 
between -8 and -1 eV. 

The obtained magnetic structure is almost antiferro- 
magnetic (without spin moment per unit cell) with nearly 
the same spin moment per Fe(l) 2+ and Fe(2) 2+ as well 
as per Fe(l) 3+ and Fe(2) 3+ cations. Using the moment 
populations in Table [IJ the calculated net moment is 
~ 0.04 \xb per Fe atom, in exact agreement with the 
experimental value fRef. Illl) . 

The charge order obtained by LSDA+L7 in 2a x b x c 
P2\/c supercell is consistent with observed enlargement 
of the (3 angle below the transition and coincides with 
c harg e ordering scheme proposed earlier by Attfield et al. 
infill It is described by the sloping 2+ and 3+ Fe cation 
lines alternately stacked along a-direction, and could be 
considered as a quasi one dimensional analog of the Ver- 



wey CO model in the pyrochlore lattice of Fe3 04. Us- 
ing the same U and J values we perform additional self- 
consistent LSDA+C/ calculations for P2\jc unit cell as 
well as for double and triple along a-direction P2\/c su- 
percells. But only self-consistent solutions with larger 
value of the total energy or with substantial magnetic mo- 
ment per unit cell, which contradicts to the experimental 
data, were found. Also we found that other charge ar- 
rangements in 2a x b x c P2±/c supercell are unstable, and 
the stable one coincides with the CO found previously. 
Thus, the CO obtained for certain value of U and J does 
not depend on the initial charge arrangemet. It is not 
possible to check all possible CO arrangements including 
more complex CO scenarios, but our results consistently 
indicate that the obtained CO solution is more favourable 
than other simple alternatives, and is the ground state of 
Fe2 0B03 in the low-temperature phase. 

Because of the small monoclinic distortion of the low 
temperature P2\/c structure the distances between the 
nearest Fe(l) cations become slightly different. It is im- 
portant to note that in the calculated CO pattern the 
shortest (2.957 A) is the distance between the pairs of 
equally charged Fe(l) cations whereas the distance be- 
tween Fe(l) 2+ and Fe(l) 3+ cations turns out to be larger 
(2.961 A). On the contrary, the pairs of Fe(l) and Fe(2) 
sites with the shortest distance between them are occu- 
pied by the cations of different valency. Thus, assuming a 
simple ionic model with the Coulomb interaction between 
Fe ions only the obtained CO pattern contradicts to the 
requirement of the minimal electrostatic energy. Indeed, 
the comparison of the Madelung energies shows that the 
energy lowers (about 0.07 eV per unit cell) if the near- 
est Fe(l) sites are occupied by Fe(l) 2+ and Fe(l) 3+ ions. 
This observation suggests that electron-lattice coupling 
rather than electrostatic repulsions drives the charge or- 
dering in Fe20B03; the same conclusion was found for 
the charge order in Fe304i2i2ii£ 

Although the corresponding total 3c? charges difference 
(0.34e) and disproportion of the total electron charges 
inside the atomic spheres of Fe 2+ and Fe 3+ cations 
(0.24e) is small, an analysis of occupation matrices of 3d 
Fe(l)/Fe(2) minority/majority spin states confirms sub- 
stantial charge separation. Thus, as shown in Tabled one 
of the ti g states of Fe(l) 2+ and Fe(2) 2+ cations is almost 
completely filled with the occupation numbers n w 0.9, 
whereas the remained two ti g orbitals of the Fe 2+ cations 
have significantly smaller population of about 0.1. Ac- 
cording to Ref. 1(J we define an order parameter as the 
largest difference between Fe 2+ and Fe 3+ ti g populations. 
While due to strong static "screening" effects^ the order 
parameter introduced as the total 3d charge difference be- 
tween 2+ and 3+ Fe cations is ill defined, the well-defined 
order parameter is the difference of tig occupancies for 
Fe 3+ and Fe 2+ cations, which amounts to 80% of ideal 
ionic CO model and clearly pronounces the existence of 
CO below the transition. The occupation matrices anal- 
ysis shows that the change of the ti g occupations caused 
by the charge ordering is very effectively screened by the 
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TABLE I: Total and /-projected charges, magnetic moments, 
and occupation of the most populated t^g orbitals calculated 
for inequivalent Fe atoms in the low-temperature P2i /c phase 
of Fe 2 OB0 3 . 



Fe ion 


1 


q s 


Qp 


Qd 


M (>b) 


t2 g orbital 


n 


Fe(l) J+ 


6.90 


0.40 


0.55 


5.95 


-4.20 




0.10 


Fe(l) 2+ 


7.12 


0.35 


0.50 


6.27 


-3.65 




0.91 


Fe(2) 3+ 


6.79 


0.38 


0.54 


5.86 


4.33 




0.09 


Fe(2) 2+ 


7.04 


0.34 


0.50 


6.21 


3.69 


dxy]' 


0.89 



TABLE II: The averaged Fe-O distances in the plane of t2 g 
orbitals for P2i/c structure of Fe20B03. d xy approximates 
to the doubly-occupied orbital of the 3d 6 Fe 2+ states. 



Fe atom 


orbital 


rforb. (A) 


d av . (A) 


Fe(l) 




2.111 


2.085 




dyz 


2.076 






d Z x 


2.069 




Fe(2) 


d X y 


2.109 


2.082 




dyz 


2.083 






d z x 


2.055 





rearrangement of the other Fe electrons. Thus, signifi- 
cant contribution to the charge screening is provided by 
Fe e g states due to relatively strong a bonds with 2p O 
states and, as a result, appreciable contribution to the 
occupied part of the valence band. 

The occupied t2 g states of Fe 2+ cations are predomi- 
nantly of d xy character in the local cubic frame (accord- 
ing to that we later mark the orbital as d xy orbital) . This 
is illustrated in Fig. O which shows the angular distribu- 
tion of the majority and minority spin 3c? electron den- 
sity of the Fe(2) and Fe(l) cations, respectively^ Thus, 
occupied Fe 2+ and unoccupied Fe 3+ cations are ordered 
alternately within the chain which is infinite along a- 
direction. The angular distribution of charge density of 
the Fe(l) and Fe(2) cations, which correspondingly be- 
longs to different Fe-ribbons being formed a cross in the 
Fe20B03 structure projected on (100) plane (see Fig.QJ 
is shown in Fig. An analysis of interatomic distances 
CTable lTH) shows that the average Fe(2)-0 distance (2.109 
A) in the plane perpendicular to one of the diagonals of 
the distorted Fe(2)06 octahedron is considerably larger 
than average distances in the other two planes (2.055 and 
2.083 A). It turns out that the occupied Fe(2) t2 g major- 
ity spin orbital is the one oriented in the plane with the 
largest average Fe(2)-0 distance. The same is also true 
for the Fe(l) ion but in this case the variation of the 
average Fe(l)-0 distances is smaller (2.111 vs 2.069 and 
2.076 A) and, as a consequence, the out-of-plane rotation 
of the occupied t2 g minority spin orbital is stronger. 




FIG. 5: The angular distribution of the majority and minor- 
ity spin 3d electron density of the Fe(2) and Fe(l) cations, 
respectively, within Fe-ribbon. The size of orbital corre- 
sponds to its occupancy. Oxygen atoms are shown by small 
spheres. X-Y-Z coordinate system corresponds to the local 
cubic frame. 




FIG. 6: The angular distribution of the majority and minor- 
ity spin 3d electron density of the Fe(2) and Fe(l) cations, 
respectively, from different Fe-ribbons,— The size of orbital 
corresponds to its occupancy. The frame of four Fe(l) atoms 
from the Fe-ribbon presented in Fig. is shown by dashed 
lines. 



V. EXCHANGE COUPLING CONSTANTS 

Using the LSDA+C/ method the exchange interac- 
tion parameters have been calculated via the variation 
of ground state energy with respect to the magnetic- 
moment rotation angle&i In Table ITTT1 we have shown the 
total set of different intraribbon exchange parameters as 
well as a contribution of different subbands into exchange 
interactions. Ji represents the effective pair exchange 
interaction between Fe atoms with effective Heisenberg 
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TABLE III: Total and partial intraribbon exchange interac- 
tion parameters are shown. The values are given in kelvin. 
The spatial representation of all these exchanges is schemati- 
cally presented in Fig. |7| 



Fe(2) 



Ji 


tig — tig 


t2g - e g 


e g e g 


Total 


Jl 


4 


-22 


-20 


-38 


■h 


-7 


-38 


1 


-44 


J 3 


3 


-19 


-16 


-32 


Ji 


4 


-72 


-1 


-70 


J s 


17 


-69 


-16 


-69 


J 6 


53 


-34 


-13 


6 


Jl 


-1 


-11 


-4 


-16 


■is 


-13 


-120 


-4 


-136 


Jo 


3 


-22 


-20 


-39 



TABLE IV: Total and partial interribbon exchange interac- 
tion parameters are shown. The values are given in kelvin. 
The spatial representation of all these exchanges is schemati- 
cally presented in Fig. [H] 



Ji 


t2g - tig 


t 2g - e g 


e g e g 


Total 


Jio 


-12 


-114 


-104 


-229 


Ju 


-16 


-140 


-28 


-183 


Jl2 


-17 


-141 


-52 


-209 


Jl3 


-3 


-52 


-2 


-56 


Jl4 


-22 


-165 


-20 


-207 


JlS 


-12 


-115 


-19 


-147 



Hamiltonian H = — J2i>j Jij e i ' e ji where e, and ej are 
magnetic moment unit vectors at site i and j. Positive 
(negative) values of J correspond to the ferromagnetic 
(antiferromagnetic) coupling between Fe sites. The spa- 
tial representation of all these exchanges is schematically 
presented in Fig. [7] Surprisingly, only the exchange inter- 
action parameter between Fe(2) 2+ and Fe(2) 3+ cations is 
ferromagnetic with relatively small value of Jq = 6 K. In 
contrast, the nearest sites in quasi-one-dimensional Fe(l) 
chain are coupled antiferromagnetically with noticeably 
larger exchange absolute value of |Js| = 69 K. Furter- 
more, the exchange parameters between the nearest sites 
of two Fe(l) chains are relatively strong and antiferro- 
magnetic (see J2, J4, and Jg in Table ITTT|l . Therefore, 
the Fe(l) sublattice is highly frustrated, while the rela- 
tively weak frustrations in the Fe(2) sublattice consider- 
ably reduce ferromagnetic interaction within Fe(2) chain. 
Also it is interesting to note that relatively strong fer- 
romagnetic intrachain interaction between ti g subbands 
of Fe(2) 2+ and Fe(2) 3+ cations (see J 6 in Table is 
strongly suppressed by the substantial antiferromagnetic 



t-2g 



e g and e g 



e g exchange. 



On the other hand, the interribbon exchange inter- 
action parameters between Fe(l) and Fe(2) atoms are 
considerably larger. The values of these interactions are 
shown in Table ITVl whereas the spatial representation is 
schematically presented in Fig. |HJ Thus, the exchange 
parameters between Fe(l) 3+ and Fe(2) 2+ cations are an- 
tiferromagnetic with values of Jio =-229 K and J14 —- 




FIG. 7: The sketch of the arrangement of exchange inter- 
action parameters within the ribbon of iron atoms. Open 
circles correspond to Fe 3+ , while Fe 2+ cations are noted by 
the closed circles. The spin moment direction on each Fe site 
is shown by an arrow. 




FIG. 8: The sketch of the arrangement of interribbon ex- 
change interaction parameters. Open circles correspond to 
Fe 3+ , while Fe 2+ cations are noted by the closed circles. The 
spin moment direction on each Fe site is shown by arrow. Note 
that Jii exchange parameters presented here have the same 
total values as Ji, while subband contributions are different. 



207 K (see Table IIV|I . Such an appreciable difference 
between Jiq and J14 arises from geometry since the for- 
mer is due to the superexchange interaction between Fe 
ions linked by shortest bonds to a common O ion. It 
seems that the geometrical reason is also responsible for 
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decreasing of absolute value of the exchange interactions 
between Fe(l) 2+ and Fe(2) 3 + cations from |J 12 | = 209 
K to |Ji 5 | = 147 K (see Table EJ. Also it is interesting 
to note that the exchange interaction between Fe(l) 3+ 
and Fe(2) 3+ cations is considerably larger than between 
Fe(l) 2+ and Fe(2) 2+ (J n and J 13 , respectively). We find 
that the interribbon exchange interactions play predom- 
inant role and determine the whole L-type ferrimagnetic 
spin structure below T c in contrast to the ferromagnetic 
intrachain order due to d 5 - d 6 superexchangeiii*i£ 

VI. SUMMARY AND CONCLUSIONS 

In summary, in the present LSDA+C/ study of the low- 
temperature P2\/c phase of Fe20B03 we found a charge 
ordered insulator with an energy gap of 0.39 eV. While 
the screening of the charge disproportion is so effective 
that the total 3d charge separation is rather small (0.34), 
the charge order is well pronounced with an order pa- 
rameter defined as a difference of t2 g occupancies of 2+ 
and 3+ Fe cations (0.8). The occupied Fe 2+ and Fe 3+ 
cations are ordered alternately within infinite along a- 
axis chains of Fe atoms. This result is remarkable in view 
of the absence of directly observed CO atomic displace- 
ments in the experimental coordinates, and demonstrates 
the utility of the LSDA+J7 method as an aide to exper- 



imental studies of CO structures. However, the charge 
order obtained by LSDA+C7 is consistent with observed 
enlargement of the p angle and coincides with charge 
ordering scheme proposed earlier by Attfield et alii It 
seems certain that Fe20B03 is charge ordered below T co , 
and the absence of the long range charge ordering from x- 
ray, neutron or electron diffraction arises from formation 
of charge order within small domains, which have been 
termed "Wigner nanocrystals" , 22 Thus, the superstruc- 
ture peaks are too weak and broad to be observed against 
background in diffraction patterns, whereas the observed 
long range monoclinic lattice distortion can arise despite 
a large concentration of defects as these preserve the di- 
rection of the monoclinic distortion, but do not propa- 
gate the coherent doubling of the lattice periodicity. An 
analysis of the exchange interaction parameters obtained 
by LSDA+J7 method inevitably results in predominance 
of the interribbon exchange interactions which determine 
the whole L-type ferrimagnetic spin structure below T c , 
in contrast to the ferromagnetic intrachain order due to 
d 5 - d 6 superexchange proposed earlier in Ref. 0. 
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